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^ \ We carry out a general-relativistic global linear stability analysis of the 

amassed carbon fuel on the surface of an accreting neutron star to determine 
CN \ the conditions under which superbursts occur. We reproduce the general obser- 
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vational characteristics of superbursts, including burst fluences, recurrence times, 



and the absence of superbursts on stars with accretion rates M < O.lMEdd; where 
^ \ MEdd denotes the Eddington limit. By comparing our results with observations, 

we are able to set constraints on neutron star parameters such as the stellar radius 
and neutrino cooling mechanism in the core. Specifically, we find that accreting 



neutron stars with ordered crusts and highly efficient neutrino emission in their 
cores (due to direct URCA or pionic reactions, for example) produce extremely 
energetic (> 10 44 ergs) superbursts which are inconsistent with observations, in 
agreement with previous investigations. Also, because of pycnonuclear burning 
of carbon, they do not have superbursts in the range of accretion rates at which 
superbursts are actually observed unless the crust is very impure. Stars with 
less efficient neutrino emission (due to modified URCA reactions, for example) 
produce bursts that agree better with observations. Stars with highly inefficient 
neutrino emission in their cores produce bursts that agree best with observations. 
Furthermore, we find that neutron stars with large radii (R ~ 16 km) produce 
very energetic superbursts that conflict with observations, even if the core neu- 
trino emission mechanism is highly inefficient. Superb urst characteristics are 
quite sensitive to several other parameters as well, most notably the composition 
of the accreted gas, concentration of carbon in the ignition region, and degree 
of crystallization of the crust. All systems that accrete primarily hydrogen and 
in which superbursts are observed show evidence of H- and He-burning delayed 
mixed bursts. We speculate that delayed mixed bursts provide sufficient amounts 
of carbon fuel for superbursts and are thus a prerequisite for having superbursts. 
We compare our global stability analysis to approximate one-zone criteria used by 
other authors and identify a particular set of approximations that give accurate 
results for most choices of parameters. 
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Subject headings: dense matter — nuclear reactions, nucleosynthesis, abundances 
— stars: neutron — X-rays: binaries — X-rays: bursts 

1. Introduction 

Superbursts are energetic thermonuclear flashes on the surfaces of accreting neutron 
stars thought to be caused by unstable carbon burning (Woosley & Taam 1976; Taam & 
Picklum 1978). They are similar to their hydrogen- and helium-burning Type I X-ray burst 
counterparts (which we refer to as "normal" bursts) in that they have thermal spectra, fast 
rises, and roughly exponential decays (see Strohmayer & Bildsten 2004 for a review). Super- 
bursts distinguish themselves, however, by their much larger fluences and longer recurrence 
times. Superbursts are a fairly new observational phenomenon: The first superburst was 
discovered in the system 4U 1735-444 by Cornelisse et al. (2000). Since then, eight more 
superbursts have been discovered in six other sources (see Kuulkers 2004 for a review). All 
nine superbursts had observed energies of ~ 10 42 ergs and were detected in systems with 
accretion rates between 10% and 30% of the Eddington limit. Recently, several more super- 
burst candidates were observed in GX 17+2, which accretes at a rate near the Eddington 
limit (in't Zand et al. 2004a,b). Superburst recurrence times are not well constrained, though 
three have been observed within 4.7 years from the system 4U 1636-536 (Wijnands 2001; 
Strohmayer & Markwardt 2002; Kuulkers et al. 2004). 

Previous theoretical studies of superbursts (Brown & Bildsten 1998; Cumming & Bild- 
sten 2001; Strohmayer & Brown 2002) have been quite successful at modeling the general 
properties of superbursts. These models demonstrate that a large carbon-enriched layer 
on the surface of an accreting neutron star can indeed produce a thermonuclear burst due 
to unstable carbon burning if the accretion rate M > O.lMEdd- The resulting energetics 
and recurrence times of the thermonuclear flashes are roughly consistent with observations. 
However, these models employ an approximate criterion to estimate the column depth at 
which the unstable carbon ignition will occur. Specifically, they define ignition to occur at 
the depth at which the condition de nuc /dT > de coo i/dT is satisfied, where e nuc is the nu- 
clear energy generation rate and e coo i is an approximation to the cooling rate. Furthermore, 
these models integrate the stellar structure and time evolution equations only down to the 
superburst ignition region, applying physically-motivated but not necessarily self-consistent 
boundary conditions at both ends of their computational domains. However, the long char- 
acteristic accretion timescale for superbursts means that the thermal diffusion depth is quite 
large, often down to the core. Thus, to accurately model superbursts, one must solve for the 
thermal and hydrostatic profiles of the entire crust, well below the ignition region. Recently, 
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Brown (2004) constructed an improved model that integrates the equations all the way to 
the core and solves for the inner temperature boundary condition using a method that is 
more self-consistent than previous studies. However, he uses the same approximate ignition 
criterion described above and assigns a fixed outer temperature boundary condition. 

Some authors have attempted to go beyond these simple approximations in modeling 
normal Type I X-ray bursts. Fushiki & Lamb (1987) were the first to discuss global per- 
turbations of the steady-state configuration of the accreted layer. Their analysis was rather 
crude, however, since it assumed a constant temperature perturbation as a function of depth 
and did not treat perturbations below the accreted layer. Narayan & Heyl (2003, hereafter 
Paper I) developed this approach much further by carrying out a full global linear stability 
analysis to determine the ignition conditions of normal Type I X-ray bursts, thereby putting 
thermonuclear burst theory on a more rigorous footing. However, their code has several 
limitations that, while relatively inconsequential for the study of normal bursts, are inade- 
quate for the application to superbursts. Explicitly, the model omits carbon burning, is not 
fully general-relativistic, and is unable to integrate the stellar structure and time evolution 
equations below the depth at which neutron drip occurs. In the present investigation, we 
remove these limitations by significantly improving and expanding the burst model of Paper 
I. We then apply the new model to superbursts. 

The motivation of this present study is to develop a mathematically more rigorous and 
self-consistent stability analysis to determine the physical conditions under which super- 
bursts occur. Our goals are to compare and contrast our results to preceding theoretical 
work, determine the effects of a wide range of neutron star parameters upon superburst 
characteristics, and put the theory of superburst ignition on a more rigorous footing to en- 
able quantitative comparisons with observations. Some of the results of this research were 
previously reported in Cooper & Narayan (2004). 

We begin the paper in §2 with a description of our numerical model, highlighting the 
additions to and improvements of the original model of Paper I. In §3, we discuss the thermal 
profile of the accreted layer and crust of the neutron star. In particular, we illustrate the 
effects on the thermal profile of varying several physical parameters such as the accretion rate, 
accreted gas composition, accreted layer composition, impurity concentration, degree of ion 
crystallization in the crust, and neutrino emission mechanism in the core. The consequences 
of these parameters upon the resulting superburst characteristics, notably the energetics 
and recurrence time, are significant, as we describe in §4. In §5 we compare our results with 
observations. We compare and contrast the results of our global linear stability analysis to 
those of the one-zone approximation in §6, and we conclude with a summary in §7. 
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2. The Model 

In this section, we outline the basic theoretical model, emphasizing the improvements 
over the original model described in Paper I. However, we will only briefly review the stability 
analysis procedure. The reader is encouraged to refer to Paper I for details. 



2.1. Governing Equations 

We assume that gas accretes spherically onto a compact object of gravitational mass M 
and areal radius R at a rate M, where M is the rest mass accreted per unit time as measured 
by an observer at infinity. We consider all physical quantities to be functions of S, which we 
define as the rest mass of the accreted gas as measured from the top of the accreted layer 
divided by AixR 2 . Near the stellar surface, S is properly interpreted as the column density. 
We use d/dt and <9/<9£ to represent the Eulerian time and spatial derivatives, respectively, 
and d/dt for the Lagrangian derivative following a parcel of gas: 

d = d_ J4_ /c ,d_ 

dt dt AnR 2 d£ ' y } 

Here <3> is the metric function (Misner et al. 1973), which reduces to the gravitational potential 
in the Newtonian limit. The term e - *^ 2 is equal to the redshift 1 + z and relates time in 
the local frame to time at infinity. We include hydrogen, helium, and carbon burning for 
the nuclear energy generation rates. Thus, we describe the composition of the gas by the 
hydrogen mass fraction X, helium fraction Y, CNO fraction Zqtsso and heavy element fraction 
Z = 1 — X — Y — -Zcno, where Z refers to all metals other than CNO. 

The stellar structure and time evolution of the accreting gas are governed by a set of 
nine partial differential equations (Thorne 1977; Brown 2000): 

* = _^l_ (1 _!^)./ 2> (2) 

oh m b nr 2 rc 2 

dm _ 4ttR 2 p ^ 2Gni y /2 
<9E rribn rc 2 ' 

"ayT = 4 1 + 2-1 r ' > W 

8P = Grr^ + f- + 4^> _ 
on rribnr 4 pc 2 mc 2 rc 2 

e- 2 * /c2 ^(^e 2 */ c2 ) = -T^ - (e H + 6 He + e c + e N - e v \ (6) 
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-*/c 2 An>*/ C ^ = 3R Ff>K (7) 
<9L l ' lQam b nW {) 

dX _ e H , 
~dT ~ (8) 

dY €h CRe /q\ 

dt ~ E* nc E* c [ } 

Note that equation (5) is the Tolman-Oppenheimer-Volkoff equation of hydrostatic equilib- 
rium. In these equations m is the interior gravitational mass, p is the mass density (such 
that pc 2 is the energy density), m h is the mass of one baryon, n is the baryon number den- 
sity, r is the areal (Schwarzschild) radius, P is the pressure, F is the energy flux, T is the 
temperature, s is the entropy per unit mass, en, eHe, £c, an d cn are the energy generation 
rates due to hydrogen burning, helium burning, carbon burning, and deep crustal heating, 
e u is the energy loss rate due to neutrino emission, k is the opacity, and E^ e , and Eq are 
the total nuclear energies released per unit mass of hydrogen, helium, and carbon burned, 
respectively. 

Note that in the post-Newtonian stellar structure equations, the partial derivatives of the 
physical parameters are taken with respect to the Eulerian variable L, while the derivatives 
should be taken with respect to a Lagrangian variable. However, we solve the equations in 
quasi-steady state. Specifically, to find the equilibrium configuration, we set d/dt — and 
solve the following set of ordinary differential equations 

*_*_(!_ ^)./» 
dL mi ) nr z rc z 

dm AnR 2 p 2Gm 



1 - ^r) 1/2 , (12) 



dL rribn v rc 2 

d$ GmR 2 , 4tit 3 P w 2Gm . 1/9 

-JV = 4 1 + 2" 1 2- ' > ( 13 ) 

dL m&nr 4 mc z rc z 

dL rribnr* pc z mc z rc z 

e-^(^) = -T^e-/^ - ( £ „ + £ „. + £c + fN - «.), (15) 

AfTp $ / c2 l - m2 Fpn 

dL 1 J 16am 6 nW 1 j 



dX _47i\R 2 $/c 2 e H 
dE ~ aF 6 



(17) 
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dY _ AnR 2 z, /c 2 ( e H _ 

dE M l £ H E*J> [ ' 

d^cNO _ 4:irR 2 g /c 2 . e He _ _ec_x , , 

Clearly, since the compact object continuously accretes matter from its companion, the 
system is never precisely in equilibrium. However, the accretion timescale is longer than that 
of most of the relevant physical processes, so the steady-state approximation is acceptable. 
In this case, E becomes a good Lagrangian variable, and so the use of E as our independent 
variable is justified. 

Stable and/or unstable hydrogen and helium burning produces most of the carbon that 
ultimately triggers a superb urst. As of this writing, each system with an observed superburst 
also undergoes normal Type I X-ray bursts (Kuulkers 2004; in't Zand et al. 2004a). The 
carbon yield resulting from both stable helium burning and unstable burning during normal 
bursts is uncertain (see §4). To account for this uncertainty, we introduce a free parameter 
Cf, the fraction of hydrogen and helium that ultimately burns to carbon. Equation (19) thus 
becomes 

dZcNO AnR 2 $ /c 2. e He e c , , . 

so at the base of the accreted layer, Zcno ~ Cf. Clearly this is an approximation. However, 
it enables us to model the composition of the accreted layer as well as possible. Specifically, 
we can model the composition of both the hydrogen- and helium-rich upper region of the 
outer crust and the carbon-rich lower region to high accuracy. Additionally, the layer in 
which the hydrogen and helium burn to carbon is most likely very narrow in the column 
density E, so this transition region is inconsequential to the thermal and hydrostatic profiles 
of the outer crust. 



2.2. Boundary Conditions 

The solution to the set of nine coupled partial differential equations (11-19) requires 
nine separate boundary conditions. Eight are applied at the photosphere (where the optical 
depth r out = 2/3) and one is applied at the crust-core interface. 

The outer boundary conditions for equations (11-13) are r = R, m = M, and $ = 
(c 2 /2) ln(l — 2GM/Rc 2 ). The outer boundary conditions for the hydrogen, helium, and CNO 
mass fraction evolution equations (17-19) are given by the composition of the accreting gas, 
such that X = X out , Y = Y ont , and Zcno = -ZcNo.out, respectively. The value for E at the 
photosphere, S out , is obtained approximately by taking the opacity at the photosphere to be 
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given by electron scattering. Thus, S out = r ou t/ (0.2(1 + X out )) gem 2 . The outer boundary 
condition for equation (14) is then given by 

GM . 2GM . U9 

= s ut^(i - -^r l/2 - (2i) 

The method we use to determine the outer boundary condition for equation (15), the 
outward flux at the stellar surface F ont , is explained in Paper I. To summarize, we assume 
a given value of F ont which, when added to the gravitational energy flux due to accretion, 
also defines a surface temperature T out . We then integrate the differential equations and 
compare the resulting temperature at the bottom, the crust-core interface, to the required 
temperature inner boundary condition. We adjust F out and repeat the integration until the 
temperature boundary condition at the bottom is satisfied to high accuracy. What makes 
our new method superior to that of Paper I is that we now integrate all the way to the 
stellar core (p ~ 2 x 10 14 gem -3 , where po denotes the rest mass density). Previously, we 
integrated only a couple of diffusion depths into the star, but we were unable to integrate 
past the neutron drip point (po ~ 4 x 10 11 gem -3 ). The long recurrence times of superbursts 
make the old method inadequate. For many calculations, the thermal diffusion depth is 
deeper than the crust-core interface, so integration to the core is necessary. 

For a given calculation, we employ one of two methods to determine the temperature 
inner boundary condition for equation (16). In our first method, we assume that we know 
the rate of neutrino emission from the core; we use two prescriptions for this, either modified 
URCA reactions (Friman & Maxwell 1979; Yakovlev & Levenfish 1995) or pionic reactions 
(Maxwell et al. 1977). Integrating to the crust-core interface gives values for the proper 
temperature T, areal radius r, interior gravitational mass to, and energy flux F. The energy 
flux, which is directed inward at the interface (F < 0), must be balanced by the neutrino 
cooling of the core. Using the neutrino luminosities L u (m,T) from Shapiro & Teukolsky 
(1983), we determine the core temperature T corc via the equation 

L u {m,T corc ) = -4vrr 2 F. (22) 

We then compare T core to T and iterate until they match. Note that we do not modify 
the formulae to account for the volume of the core. Since we do not know the equation 
of state of the core, we cannot solve for the core's proper volume. However, since the 
formulae are estimates, and since the luminosity is a very strong function of the temperature 
(L™ URCA oc T 8 and oc T 6 ), this approximation is innocuous. In our second method, 
we simply assign a value for the core temperature T core and iterate until T = T corc at the 
crust-core interface. 
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2.3. Auxiliary Conditions 

2.3.1. Equation of State 

We assume photons, electrons, ions, and free neutrons supply the pressure. To calculate 
the photon and electron contributions to the pressure, we follow the same procedure as in 
Paper I, while for the ions and free neutrons, we follow the prescriptions of Brown (2000). 
The phase of the ions, whether solid or liquid, is determined by the dimensionless coupling 
parameter 

where nj is the ion number density. In the outermost layers of the crust, where T < 1, we use 
the free energy fits of Chabrier & Potekhin (1998) for the ion equation of state. For T > 1 
we use the analytical fits of Farouki & Hamaguchi (1993) for the Helmholtz free energy to 
calculate the ion pressure. We assume the ions are in the liquid phase when 1 < T < 173 
and the solid phase when T > 173. The analytical fits are valid for a one-component plasma. 
However, in the accreted layer more than one species is usually present at a given density. 
In this case we approximate the multi-component mixture as a one-component plasma by 
making the substitutions Z — > (Z) and rij — > n/fi nuc , where fi nuc is the mean molecular 
weight per nucleus. Below the accreted layer, we assume only one species is present at 
a given density. We calculate the free neutron pressure from the compressible liquid-drop 
nuclear model of Mackie & Baym (1977). For p > 1-3 x 10 13 gcm~ 3 , we use the equation 
of state from Negele & Vautherin (1973). 

The photon entropy formula is the same as that in Paper I. We calculate the Coulomb 
portion of the ion entropy from the free energy fits of Chabrier (1993) and the ideal portion 
from the fits of Farouki & Hamaguchi (1993). The treatment of the electron and free neutron 
entropies is more difficult. The compressional heating terms in equation (15) are unlikely 
to affect the thermal profile calculation significantly (Brown & Bildsten 1998) and there- 
fore are often omitted in other investigations. However, we attempt to make a reasonable 
approximation to the compressional heating terms relevant to our work. We consider three 
entropy "regimes": ideal, nonrelativistic degenerate, and extremely relativistic degenerate. 
Define 6 = m p k B T / \h 2 n 2 ^) and x = (^/m p c)(37r 2 n) 1 / 3 , the relativity parameter (Shapiro 
& Teukolsky 1983), where m p is the particle mass, k B is Boltzmann's constant, and n is the 
particle number density. Then the expressions for the entropy of a species in each of the 
three regimes are (Landau & Lifshitz 1969) 
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Note that ^> 1 for an ideal gas and 6 < 1 for a degenerate gas. Similarly, x 3> 1 for 
a relativistic gas and x <C 1 for a nonrelativistic gas. The question is, then, how does one 
decide where to make the transitions between the three regimes? It is evident from equation 
(15) that only the entropy derivatives, not the entropy values, are important. Therefore, we 
choose the transition values of 6 and x in such as way as to make the partial derivatives of 
the entropies with respect to the variable ori continuous functions of that variable. First, 
we decide whether the species is ideal or degenerate. Thus, for a quantity q, we first find the 
transition value Otrans such that dsi^i/dq = <9snrd/<9<? at 9 = O tr ans- For all quantities q, 
©trans = 1.455. If the species is degenerate, we then decide whether it is nonrelativistic or 
extremely relativistic. We then repeat the same procedure with 6 — > x. x tTSulls = 0.5, 1, or 
2, depending on the specific quantity q. 



The hydrogen and helium nuclear energy generation rates are identical to those in Paper 
I with the exceptions that the 13 N(p, 7) 14 and 14 N(p, 7) 15 rates for e B are updated to those 
of Caughlan & Fowler (1988). 

Carbon burns through several different energetically possible reaction channels (Clayton 
1983). For simplicity, we assume that the dominant reaction is 12 C( 12 C,7) 24 Mg, so = 
5.6 x 10 17 ergg _1 (Cumming & Bildsten 2001). For ec, we use both the reaction rates and 
electron screening enhancement factors of Kitamura (2000). We interpolate the logarithms 
of the fluid and solid phase reaction rates in the same way that Kitamura (2000) interpolates 
the thermonuclear and pycnonuclear enhancement factors. This makes the energy generation 
rate a smooth function of both density and temperature. See Figure 1 for a plot of ec 
as a function of po- Note that we use the approximate criterion of Kitamura (2000) for 
the solidification of an electron-screened binary ionic mixture. Previous authors (Brown & 
Bildsten 1998; Cumming & Bildsten 2001; Strohmayer & Brown 2002; Brown 2004) have 
used the reaction rate of Caughlan & Fowler (1988) with the enhancement factor of Ogata 
et al. (1993). This energy generation rate is valid only in the liquid phase, where the 
burning is considered strictly thermonuclear. In contrast, our energy generation rate is valid 
in both the thermonuclear regime, where the crust is fluid, and the pycnonuclear regime, 
where the crust is solid. The two energy generation rates agree to sufficient precision in 



2.3.2. Nuclear Energy Generation Rates 
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the thermonuclear regime. In the pycnonuclear regime, the nuclear energy generation rate 
is essentially temperature-independent. Including this regime in the model can sometimes 
be important, especially in neutron stars with efficient core cooling mechanisms or heavy- 
element oceans. 



2.4. Deep Crustal Heating and Neutrino Emission 

For 6 x 10 11 gem" 3 < p < 3 x 10 13 gem -3 , electron captures, neutron emissions, and 
pycnonuclear reactions release ~ 1 MeV per baryon over a timescale inversely proportional 
to the accretion rate (Haensel & Zdunik 1990a,b; Brown et al. 1998). For e^, the energy 
generation rate due to this deep crustal heating, we use the formula of Brown (2000). 

We include pair, photo-, plasma, and bremsstrahlung neutrino emission processes for 
the neutrino energy loss rate e v . For the pair, photo-, and plasma processes we use the 
analytical fits of Itoh et al. (1996), and for the bremsstrahlung emission rate we use the work 
of Haensel et al. (1996) when the ions are in the liquid phase and Yakovlev & Kaminker 
(1996) when the ions are in the solid phase. 

3. Thermal Structure of the Accreted Layer and Crust 

Superburst characteristics depend sensitively upon the thermal structure of the neutron 
star crust. In this section, we identify some of the key parameters that affect the thermal 
profile. 

3.1. Accretion Rate 

In its journey from the binary companion to the neutron star core, a parcel of matter 
releases a tremendous amount of energy, and the rate of this energy release is directly pro- 
portional to the rate at which matter falls onto the stellar surface. Thus, it is not surprising 
that the thermal profile of an accreting neutron star is very sensitive to the accretion rate. To 
varying degrees of importance, the accretion rate affects the thermal profile in four different 
ways. 

Most of the energy released by an infalling parcel of accreted material is from the grav- 
itational energy released when the matter impacts the stellar surface. Most of this energy 
is radiated outward. Nevertheless, it determines the temperature at the stellar surface and 
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thereby sets a boundary condition. The calculation is rather insensitive to this temperature, 
however, for the thermal profile near the surface approaches a radiative-zero solution. As 
accretion continues, nuclear fuel accumulates and eventually burns either stably or unstably. 
The time-averaged rate of nuclear energy generation is proportional to the accretion rate. 
Although the nuclear energy per gram of accreted material released via fusion is roughly 
forty times less than that released from gravitational energy, the nuclear energy is generated 
well below the stellar surface. Thus the nuclear energy generation can have a significant 
effect upon the thermal profile of the crust, in particular the superb urst ignition region (see 
also §3.2). Additionally, continuous accretion causes both compressional heating throughout 
the crust and deep crustal heating via non-equilibrium reactions. The compressional heat- 
ing is rather small compared to other sources, and it is therefore often neglected in other 
studies. The deep crustal heating is roughly five times less than that from hydrogen and 
helium burning, but it can have a non- negligible effect upon the thermal profile of the crust, 
especially if the conductive opacity of the inner crust is large. 

We plot the temperature and flux profiles for two neutron stars accreting at different 
rates in Figure 2. The energy flux is normalized by the maximum nuclear burning energy 
flux available in the accreting gas: 

M 2GM 

^ nuc = AnR 2 ^ Rc^~^ 1 ^ 2 [X utEft + (X out + Y ont )E^ c + (X out + Y ont + Z C NO,out)-Ec]- (27) 

The parameter / acc = M/MEdd is the accretion rate normalized to the Eddington limit, where 
^Edd = 47rGM(l + z)/czn es , with n es = 0.4 cm 2 g~ 1 . The temperature outer boundary 
condition, which is shown at the left end of the left panel, is determined by the rate of 
gravitational energy liberated at the surface. Changes in the slope of the temperature profile 
are associated with localized energy sources, and they are reflected by rapid changes in 
the flux. Thus, hydrogen and helium burning causes the peak in the thermal profile and 
the large change of flux at a column depth E « 10 8 gcm~ 2 . A small amount of carbon 
burning occurs at £ ~ 10 12 -10 13 gcm~ 2 . Deep crustal heating occurs for column depths 10 15 
gcm~ 2 < £ < 10 17 gem" 2 . The flux profile shows that most of the energy generated by 
deep crustal heating is directed inward, in agreement with Brown (2000). 



3.2. Energy Generated from Hydrogen and Helium Burning 

The burning of hydrogen and helium near the surface of an accreting neutron star re- 
leases a substantial amount of energy within the star. Consequently, the thermal profile of 
the outer crust, including the superb urst ignition region, is rather sensitive to the magni- 
tude and physical location of this energy generation. To account for this, previous authors 
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set the temperature at a given column depth to coincide with estimates from investigations 
of hydrogen and helium ignition. However, the thermal profile in this region is a sensitive 
function of many variables, including the mass accretion rate (e.g., compare the two models 
in Fig. 2), stellar radius, and composition of the accreted gas. Since we include both hydro- 
gen and helium energy generation rates in our energy conservation equation, we make no 
assumptions regarding the temperature at a given depth in the accreted layer. Thus, we are 
able to determine the thermal profile of the outer crust self-consistently. Not only does this 
improve the accuracy of our calculation, but it also gives us the freedom to vary physical 
parameters such as the gas composition and stellar radius self-consistently. 

As noted earlier, all of the systems in which astronomers have observed superbursts 
exhibit normal Type I X-ray bursts as well. To do a rigorous calculation of the thermal 
profile of the outer crust, one would need to conduct a fully time-dependent calculation of 
many successive normal bursts, which is beyond the scope of this study. Our calculation 
is quasistatic, so the composition and thermal profile of the crust in the the normal burst 
ignition region is essentially computed via stable (though rapid) hydrogen and helium burn- 
ing. Since the timescale over which normal bursts occur (hours to days) is much shorter 
than the timescale over which superbursts occur (years to possibly decades), any effects that 
hydrogen and helium burning have on the thermal profile of the superburst ignition region 
will be due to the time-averaged hydrogen and helium nuclear energy generation rate. The 
time-averaged energy generation rate is the same regardless of the manner in which the fuel 
is burned. Therefore, our method should be sufficiently accurate for our purposes. 

To demonstrate the importance of hydrogen and helium burning on the thermal profile 
of the outer crust, we plot in Figure 3 the temperature as a function of column density for 
two systems with different accreted gas compositions. Hydrogen burning releases much more 
energy per gram of accreted fuel than helium burning. Therefore, the maximum temperature 
achieved in the outer crust is in general proportional to the mass fraction of hydrogen of the 
accreted gas. This can have a significant effect on superburst characteristics (see §4.1). The 
inner crust (E > 10 15 gem -2 ), however, is rather insulated from the hydrogen and helium 
burning region, so the thermal profile of the inner crust is rather insensitive to the hydrogen 
and helium burning near the surface. 



3.3. Ash Composition from Hydrogen and Helium Burning 

The heavy element composition of the accreted layer where superbursts are triggered 
is uncertain. Especially for mixed hydrogen/helium accretors, the primary cause of this 
uncertainty is due to the rp-process (Wallace & Woosley 1981), whose ashes are probably a 
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mix of elements beyond the iron peak (Schatz et al. 2001). Type I X-ray burst models of 
Woosley et al. (2004) produce nuclei with an average atomic weight (A) pa 64. Therefore, we 
choose two representative metals as our heavy elements: (i) 26 Fe, and (ii) following Cumming 
& Bildsten (2001), ^Ru. These choices likely bracket the true average atomic weight and 
electric charge of the ashes. We assume that the nuclear composition below the accreted layer 
is that of Haensel & Zdunik (1990a,b) for both cases. Haensel & Zdunik (2003) later studied 
the evolution of heavy ({A) pa 100) rp-process ashes. However, Schatz et al. (2003a) showed 
that the high temperatures reached during a superburst might induce photodisintegration 
reactions in the heavy ashes, converting them to iron group elements. Thus, we presume 
that the outer crust below the accreted layer consists chiefly of iron group elements, with 
possible impurities. 

We find that, contrary to the study of Cumming & Bildsten (2001), the composition 
of the ashes from hydrogen and helium has a negligible effect on the thermal profile of the 
crust. This is consistent with the results of Brown (2004). See Figure 4 for a plot of the 
thermal profile for neutron star crusts with different ash compositions. 



3.4. Core Temperature 

The thermal profile of both the inner crust and the outer crust below the hydro- 
gen/helium burning region is very sensitive to the temperature at the crust-core interface. 
See Figure 5 for a plot of the thermal profiles of neutron stars with different core temper- 
atures. In particular, note the slopes of the profiles beyond the hydrogen/helium burning 
region, at E pa 10 s gem" 2 . Except for neutron stars with very hot cores (T > 4 x 10 8 K), the 
temperature gradient dT/dT, in this region is negative (i.e. dT/dr is positive), which implies 
that the net energy flux is negative. This is clearly illustrated in the energy flux profiles 
in Figure 5. Previous studies of thermonuclear burst ignition have often used the result of 
Brown (2000) for the inner flux boundary condition. He found that approximately 10% of the 
energy generated in the crust through deep crustal heating flows outward. In Figure 5 this 
corresponds to F/F QUC pa 0.02. For comparison, the energy flux in Figure 5 just below the 
column depth at which the superburst is triggered ranges from F/F nnc pa +0.05 at E pa 10 11 ' 7 
gem -2 for the hottest thermal profile to F/F nnc pa —0.06 at E pa 10 14 ' gcm~ 2 for the coldest. 
Brown's results, which are accurate for the particular physical scenario he studied and which 
we are able to reproduce, are clearly not applicable for all scenarios. As we have shown, the 
thermal profile of an accreting neutron star is sensitive to many parameters, including the 
accretion rate, composition of the accreted gas, and core temperature. Therefore, one must 
use caution when implementing flux or temperature boundary conditions. In our opinion, it 
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is preferable to do self-consistent calculations, as in this work. 

In recent work, Brown (2004) has used an approximate "outer" boundary condition on 
the temperature: T = 2.5 x 10 8 K at E = 10 9 gem -2 . A quick look at Figures 2, 3, and 
5 here shows that the temperature varies considerably at this depth, depending on various 
parameters. Also, T is typically greater than 2.5 x 10 8 K at this depth. 



Electron- impurity scattering (Itoh & Kohyama 1993) and electron-phonon scattering 
(Baiko & Yakovlev 1995) determine the thermal conductivity of the solidified inner crust of 
an accreting neutron star. Ion impurities embedded in a crystalline lattice cause deviations 
in an otherwise periodic potential. Electron scattering off of these deviations reduces the 
thermal conductivity of the lattice and consequently raises the temperature of the inner 
crust. Unlike electron-phonon scattering, which is due to thermal oscillations of ions and 
is therefore temperature-dependent, electron-impurity scattering depends primarily on the 
lattice composition and thus the corresponding collision frequency is essentially temperature- 
independent. Therefore, impurities can have a significant effect on the thermal conductivity 
of the crust, especially at low temperatures (Itoh & Kohyama 1993; Ashcroft & Mermin 



The contribution of impurities to the thermal conductivity of the crust is parametrized 
by the impurity parameter (Itoh & Kohyama 1993; Brown 2000) 



where rij and Zj are the number density and charge of the jth ion species, nj = rij 
is the total ion number density, and (Z) = nj 1 J2j n j^j * s the m ean charge. Sources of 
impurities in the crust include electron captures and pycnonuclear reactions (Haensel & 
Zdunik 1990a,b), but the main source of impurities is most likely the mixture of superburst 
and/or rp-process ashes at the top of the substrate. Early studies of the composition of rp- 
process ashes implied that impurity scattering should be significant, with Q ~ 100 (Schatz 
et al. 1999). Subsequent calculations by Woosley et al. (2004) and Koike et al. (2004) showed 
that rp-process ashes consist chiefly of iron-peak nuclei, resulting in a lower Q. Schatz et al. 
(2003a) found that Q as 5.2 after a superburst is triggered. In our model, we calculate Q 
self-consistently within the accreted layer, and we adopt Q = 5.2 as our fiducial impurity 
parameter value in the substrate below the accreted layer. Figure 6 illustrates the effects of 
electron-impurity scattering on the thermal profile of an accreting neutron star. 



3.5. Ion Impurities in the Crust 



1976). 




(28) 
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3.6. Ion Crystallization in the Crust 

When calculating the conductive opacity in the inner crust of the neutron star, we 
usually assume that the ions form an ordered crystal lattice when the dimensionless coupling 
parameter T > 173 (eq. 23). In this case, the conductive opacity is usually dominated 
by electron-phonon scattering due to oscillations of the ions in the lattice (Potekhin et 
al. 1999). However, previous studies of the nuclear structure of the inner crust suggest 
that some fraction of the inner crust may in fact be disordered (Magierski & Heenen 2002; 
Magierski & Bulgac 2004). Following Brown (2004), we investigate the thermal profile of 
a completely disordered neutron star crust. We use the thermal conductivity expression of 
Itoh & Kohyama (1993) and set the structure factor (S) to unity. Additionally, we set the 
squared impurity charge ((AZ) 2 ) = (Z) 2 . This essentially sets a lower limit on the thermal 
conductivity of the crust (E. Brown, private communication). The thermal profiles of the 
crusts for neutron stars with cores that emit neutrinos via either modified URCA reactions 
or pionic reactions are shown in Figure 7. In contrast to neutron stars with crystalline crusts 
(Fig. 5), the thermal profiles now are quite insensitive to the nature of the core neutrino 
cooling mechanism. 

4. Results 

4.1. Composition of the Accreted Gas 

For systems that have exhibited a superburst and for which the composition of the 
accreted gas can be reasonably estimated, all but one accrete a mixture of hydrogen and 
helium, with hydrogen being the most abundant species by mass (Kuulkers 2004). The 
system 4U 1820-30 (Strohmayer 2000; Strohmayer & Brown 2002) is the exception. Several 
studies (Fedorova & Ergma 1989; Podsiadlowski et al. 2002; Cumming 2003) imply that 
the compact star in this system accretes a helium-rich mixture with a small hydrogen mass 
fraction X ~ 0.1. Observationally, the superburst from 4U 1820-30 is distinct, with a larger 
fluence, luminosity, and peak temperature than every other superburst observed thus far 
(Kuulkers et al. 2002b; Kuulkers 2004). Therefore, to determine the effects of accreted 
gas composition on superbursts, we choose two different elemental abundances: (i) "mixed 
hydrogen/helium", for which the mass fractions of the accreted gas are X = 0.7, Y = 0.28, 
and Z CNO = 0.016, and (ii) "helium", for which X = 0.1, Y = 0.88, and Z CNO = 0.016. 

Figure 8 shows the superburst energies and recurrence times as a function of accretion 
rate, as calculated by our model. All other parameters being equal, helium accretors require 
a larger column density of accreted gas before a superburst is triggered. Therefore, their 
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superbursts are more energetic and have longer recurrence times, in agreement with the 
observations of 4U 1820-30. Helium burning releases less energy per gram than hydrogen 
burning by approximately one order of magnitude. Consequently, for a given column density, 
the temperature at the base of the layer is lower (see Fig. 3), so a larger column of fuel must 
accumulate before a superburst can occur. This disparity is greater at high accretion rates, 
for which a lower column density is required for an instability. Less matter exists between 
the hydrogen/helium burning region and the base of the layer, so the base is less insulated 
from the burning region and therefore more sensitive to the energy generated there. 

As we have shown, the composition of the accreted gas affects superburst characteristics 
by its effect upon the thermal profile of the outer crust where superbursts are triggered. 
However, the composition may have an even greater influence on superburst characteristics 
through the carbon yield resulting from both the stable and unstable burning of the gas. 
Unfortunately, we cannot investigate this aspect of the problem with our model since we 
do not solve for the carbon yield self-consistently, but set it through the parameter Cf. We 
assume the fiducial value Cf = 0.3 for all calculations unless notified otherwise. 



4.2. Composition of the Ashes in the Accreted Layer 

The composition of the accreted layer where superbursts are triggered is uncertain. 
In particular, the mass fraction of carbon produced via stable and unstable hydrogen and 
helium burning is unknown, but it most likely has to be > 10% for mixed hydrogen/helium 
accretors (Cumming & Bildsten 2001) and > 30% for helium accretors (Strohmayer & Brown 
2002). These authors find that, for M < 0.3MEdd ; a smaller mass fraction of carbon than 
these limits will not produce a superburst. Especially for mixed hydrogen/helium accretors, 
the primary cause of the uncertainty in the carbon fraction is due to the rp-process (Wallace 
& Woosley 1981). Several research groups have studied the final products of the rp-process 
from both stable and unstable burning (Schatz et al. 1999, 2003b; Koike et al. 2004; Woosley 
et al. 2004; Fisker et al. 2004). They find that the carbon mass fraction of the ashes is notably 
below 10% for the range of accretion rates at which superbursts have been observed, which 
is problematic since such a low fraction is insufficient to trigger a superburst. The exception 
to this statement may be the system GX 17+2, in which extremely energetic Type I X-ray 
bursts were observed at accretion rates in the neighborhood of the Eddington limit (Kuulkers 
et al. 2002a; in't Zand et al. 2004a). Despite these results, we cannot rule out the possibility 
that some other process (e.g., delayed mixed bursts, see §5) may still be able to produce 
enough carbon. Due to this uncertainty, we take as a free parameter Cf, the fraction of 
hydrogen and helium that ultimately burns to carbon. 
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Figure 9 shows superburst energies and recurrence times as a function of accretion rate 
for three choices of Cf. For a given accretion rate, the superburst recurrence time will usually 
be shorter for a larger mass fraction of carbon at the base of the accreted layer because an 
instability will occur sooner if more carbon is present. To first order, the burst energy 
-EWst oc Si aycr Cf, where Si ayer is the column depth of the accreted layer when a superburst 
is triggered. However, Si ayer itself is roughly proportional to the recurrence time, so it is 
difficult to derive a general relationship between Cf and E hurst (compare the burst energies 
for C f = 0.3 and C f = 0.5 in Figure 9). In general, the value of C f does not substantially 
affect superburst energetics and recurrence times for a given accretion rate. However, the 
value of Cf does significantly affect the lower limit of the range of accretion rates at which 
superbursts occur. We find that the lower limit is roughly inversely proportional to Cf, in 
agreement with Cumming & Bildsten (2001). 

As noted in §3.3, the composition of the heavy elements in the accreted layer has a 
negligible affect on the thermal profile of the superburst ignition region. Consequently, 
for accretion rates at which superbursts are triggered, superburst energies and recurrence 
times are quite insensitive to the nuclear composition of the heavy elements in the accreted 
layer. This agrees well with the results of Brown (2004). At lower accretion rates, however, 
the greater charge of the ^Ru ions may cause the base of the accreted layer to solidify 
before a thermonuclear instability can occur. In this case, the carbon fuel burns stably via 
pycnonuclear reactions (see §2.3.2), so a superburst does not occur. Therefore, though the 
composition of the heavy elements in the ocean is unimportant with regard to superburst 
characteristics, it is important with regard to the presence or absence of superbursts at a 
given accretion rate. This dependence has not been noted in previous studies. See Figure 10 
for a plot of superburst energies and recurrence times for neutron star crusts with different 
heavy element compositions. 



4.3. Neutrino Cooling Mechanism in the Core 

The composition of dense matter in the inner cores of neutron stars is essentially un- 
known. However, this composition significantly affects the neutrino emission there. Thus, 
knowledge of which neutrino processes occur in the core can help constrain the types of 
matter that exist at such high densities. 

When we solve for the core temperature by balancing the flux flowing inward with the 
neutrino cooling, we find T core m 3 x 10 8 K for nonsuperfluid cores that emit neutrinos via 
modified URCA reactions and T corc as 2 x 10 7 K for cores that emit neutrinos via pionic 
reactions. In general, T corc oc M 1 / 8 for modified URCA cooling and T core oc M 1//6 for pionic 
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cooling. The thermal profile of much of the crust, and in particular the region in which 
superbursts are triggered, is quite sensitive to the core temperature (see Fig. 5). Thus, 
compared to stars with inefficient cooling mechanisms in their cores, neutron stars with 
efficient cooling mechanisms in their cores need to accrete more fuel in order for the carbon 
at the base of the accreted layer to ignite. Consequently, these superbursts will be more 
energetic and will have longer recurrence times, as shown in Figure 11. 

If neutrons and protons exist in the stellar core, they probably become superfluid at high 
densities (Baym et al. 1969). Baryon superfluidity drastically suppresses neutrino emission 
from modified and direct URCA processes and therefore raises the temperature in the core 
(Yakovlev et al. 1999, 2001; Yakovlev & Pethick 2004). Due to the large uncertainties in the 
composition and physical state of the matter, we do not attempt to model the superfluid 
core in detail. Rather, to study the effects of superfluidity on superburst characteristics, 
we simply fix the temperature at the crust-core interface at several values and carry out 
the calculations. Figure 12 shows the results. Neutron stars with hotter cores exhibit 
superbursts with smaller fluences and shorter recurrence times. In addition, we see that for 
these temperatures the lower limit of the accretion rate range over which superbursts occur 
increases as the core temperature increases. The high crust temperature causes the carbon 
to burn stably before an instability is triggered, in agreement with Cumming & Bildsten 
(2001). 

Impurities in the crust of a neutron star can significantly affect the thermal profile of 
the superburst ignition region, especially in stars with cold cores (see §3.5). Consequently, 
superburst energies and recurrence times are quite sensitive to the concentration of impuri- 
ties in the crust. Figures 13 and 14 show some results. We see that there is little difference 
between Q = and Q = 5.2 when the core radiates neutrinos via modified URCA reactions, 
but Q = 100 induces a substantial change in the results. However, recent investigations sug- 
gest that impurity scattering in the crusts of accreting neutron stars that exhibit superbursts 
is rather insignificant (Schatz et al. 2003a). If the impurity concentration Q ~ 1 in the crust 
of a neutron star with a core that emits neutrinos via pionic or direct URCA reactions, then 
the carbon fuel will solidify and burn stably via pycnonuclear reactions at lower accretion 
rates. 

We have shown in §3.6 that, for a completely disordered neutron star crust, the thermal 
profile (including that of the superburst ignition region) is highly insensitive to the core 
neutrino emission mechanism. In this case, superburst energies and recurrence times should 
be insensitive to the core cooling mechanism, as confirmed in Figure 15. Thus we verify the 
result of Brown (2004) that superburst energies and recurrence times from neutron stars with 
highly efficient core neutrino cooling mechanisms may still be consistent with observations 
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if the crust is disordered. However, in this case we find that superbursts should occur 
even at relatively low accretion rates M < O.lMEdd, whereas observations indicate a cutoff 
at M £3 O.lMEdd- Normally, the carbon would burn stably at these low accretion rates 
because the energy generated would be efficiently transported away from the burning region. 
However, the low thermal conductivity due to the disordered lattice inhibits the flow of energy 
away from the carbon-burning region. Therefore, even a low carbon energy generation rate 
can initiate a thermonuclear instability. Note that in deriving these results we assume that 
the entire crust is completely disordered, which is clearly an extreme situation. Further 
investigations into the nuclear structure of neutron star crusts are necessary to determine 
the significance of these results. 

4.4. Stellar Radius 

The radius of a neutron star depends quite sensitively on the core equation of state, 
but it is virtually independent of the stellar mass. Accurate measurements of the radius 
to within about one kilometer can potentially constrain the equation of state (Lattimer & 
Prakash 2001). To demonstrate the effects of the stellar radius on superburst characteristics, 
we choose three different values, R = 16.4, 10.4, and 6.5 km, which likely bracket the true 
radii of neutron stars. We find that stars with larger radii have more energetic superbursts 
and longer recurrence times at a given accretion rate. They also have superbursts at lower 
accretion rates. See Figure 16. 

At a given accretion rate, a neutron star with a larger radius requires a larger column 
density of fuel in order for a superburst to be triggered. This is a result of the lower 
gravitational acceleration near the stellar surface. Thus the effect of radius on superburst 
recurrence times is twofold. Not only is the accretion rate per unit area smaller for a larger 
star, but also the amount of fuel per unit area required for a superburst to occur is larger. The 
effect of radius on superburst energetics is even stronger. In addition to the factors stated 
above, the total surface area of the star is larger, so the total amount of fuel available to burn 
at a given column depth is greater. Figure 16 thus suggests that superburst observations 
may be useful for constraining neutron star radii. 

5. Comparison with Observations 

As noted in §1, the nine superbursts that have been observed, excluding those from 
GX 17+2, have integrated photon fluxes of 10 42 ergs. The paucity of data makes the 
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recurrence time of superbursts difficult to determine, but observations imply a recurrence 
time of ~ 1-2 years. All nine superbursts occurred in systems with accretion rates between 
10% and 30% of the Eddington limit. Additionally, several superburst candidates have been 
observed in the near-Eddington accretor GX 17+2 with energies of ~ 5 x 10 41 ergs. A 
successful theoretical model of superbursts must explain these facts. 

We begin by discussing the M range over which superbursts are observed. The lack 
of superbursts for accretion rates M < O.lMEdd arises naturally in our model and is rather 
simple to explain theoretically. As discussed in §4.3, at low M, either the carbon fuel 
burns stably via thermonuclear reactions before an instability is triggered (for high core 
temperatures, Cumming & Bildsten 2001) or the crust solidifies in the superburst ignition 
region and the carbon burns stably via pycnonuclear reactions (for low core temperatures). 
In either case, there is a cutoff of superbursts at low values of M. 

The above discussion is invalid if the crust is highly disordered, since then the thermal 
conductivity is low, causing the region in which superbursts are triggered to be hotter than 
if the crust were ordered. Also, when the carbon ignites, the low thermal conductivity 
inhibits the diffusion of the nuclear energy generated there. The nuclear energy generation 
rate exceeds the rate at which thermal conduction can cool the region, so a thermonuclear 
runaway ensues and a superburst is triggered. Therefore, as seen in Figure 15, superbursts 
occur down to accretion rates below O.lMEdd- The lack of observed superbursts for accretion 
rates M < O.lMEdd then needs to be explained. One possibility is that, for some reason, e.g., 
lack of delayed mixed bursts (see below), carbon is not produced in sufficient quantities to 
fuel a superburst. Another explanation could be that at low accretion rates the recurrence 
times are so long that astronomers simply have not observed any of the systems long enough 
to see superbursts. 

The absence of superbursts for 0.3MEdd ^ M < l.OMEdd is not well understood. Every 
theoretical study, including our own, suggests that superbursts should occur more easily at 
higher mass accretion rates, assuming a sufficient amount of carbon is present. The peak 
luminosities of superbursts are often below the Eddington limit. Therefore, one possible 
explanation is that none has been observed in this range of accretion rates because there is 
little contrast between the peak burst luminosity and the accretion luminosity (Strohmayer 
& Bildsten 2004). On the other hand, the recurrence time is inversely proportional to the 
accretion rate to first order, so one would expect to observe more superbursts at higher 
accretion rates. 

Another possible explanation involves the rp-process. Perhaps the rp-process leaves 
behind enough unburnt carbon to produce superbursts only for M < 0.3MEdd, and it burns 
too much carbon at higher mass accretion rates. Both the critical amount of carbon needed 
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for superbursts and the amount of unburnt carbon left by the rp-process decrease with 
increasing M. Perhaps the two dependencies conspire to permit superbursts only for M < 
0.3MEdd- One problem with this explanation is that all studies of the rp-process carried out 
so far predict too little carbon to produce superbursts at all M. 

The above proposals do not explain the observed superbursts from GX 17+2, which 
accretes at M ~ M^dd- If one can observe superbursts in a system accreting at nearly the 
Eddington limit, than presumably one could observe superbursts in systems accreting at any 
lower rate as well. Furthermore, as discussed in §4.2, the minimum mass fraction of carbon 
needed to trigger a superburst is inversely proportional to the accretion rate. However, the 
amount of carbon left behind by stable burning via the rp-process falls even more rapidly 
with increasing M. 

Yet another possible explanation involves ideas described in Paper I. In that paper, 
Narayan and Heyl showed that at high accretion rates (M > O.lMEdd) normal Type I X- 
ray bursts occur in a unique regime that they refer to as "delayed mixed bursts." In these 
systems, a large fraction of the hydrogen and helium fuel burns stably to carbon before 
the instability is triggered. This stable burning explains the observations of van Paradijs 
et al. (1988) who found that, for systems that accrete at a high rate and exhibit normal 
bursts, the quantity a, the ratio of the energy released between bursts to the energy released 
during a burst, rises dramatically. Such an increase of a is nicely reproduced by the model 
described in Paper I. Furthermore, in't Zand et al. (2003) showed that systems that exhibit 
superbursts generally have a values that are significantly greater than systems that do not 
exhibit superbursts. 

Delayed mixed bursts may be the source of the substantial amount of carbon needed to 
trigger a superburst. As noted earlier, every system in which a superburst has been observed 
exhibits normal bursts as well. We tentatively suggest that the occurrence of delayed mixed 
bursts is a necessary condition for superbursts in systems for which the accreted material 
is predominantly hydrogen. The stably burned material forms a thick layer of carbon, 
and the delayed burst that finally occurs burns above this layer, leaving the carbon largely 
unaffected. This would explain the absence of observed superbursts in systems with accretion 
rates 0.3M Ed d < M < 1.0M Ed d- These systems do not exhibit normal Type I bursts, so 
the hydrogen and helium fuel burns stably via the rp-process. According to estimates in 
the literature, the rp-process leaves behind too little carbon for a superburst. Note that 
the existence of normal bursts in GX 17+2 is not understood theoretically if the accreted 
material is indeed predominantly hydrogen. However, Kuulkers et al. (2002b) find that 
a > 1000 in this system, which implies that the normal bursts are delayed bursts. We 
emphasize that this explanation is merely a hypothesis, and it is not without issues. In 
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particular, we do not know if the carbon produced during the stable burning phase of a 
delayed burst system would survive when the burst is triggered above it. More research is 
needed on this interesting problem. 

Delayed mixed bursts may also explain the absence of superbursts for M < 0.1M E dd 
in models with disordered crusts. These models predict superbursts down to relatively 
low accretion rates, provided they have enough carbon. However, Paper I showed that at 
accretion rates below O.lMEdd the normal Type I bursts are prompt bursts. Such bursts 
are expected to burn all the fuel to elements much heavier than carbon, and so their ashes 
should not support superbursts. 

We next consider the energies and recurrence times of superbursts. Leaving aside models 
with disordered crusts, Figure 11 shows that neutron stars with ordered crusts and highly 
efficient neutrino emission mechanisms in their cores have much more powerful superbursts 
and longer recurrence times than are observed. Furthermore, if the impurity parameter 
Q ~ 1, these stars do not exhibit superbursts at accretion rates M < 0.3M E dd, where all 
superbursts except those from GX 17+2 have been observed (see Fig. 14). Therefore, we 
conclude that accreting neutron stars with highly efficient neutrino emission in their cores 
(due to direct URCA or pionic reactions, for example) are inconsistent with superburst 
observations, in agreement with Brown (2004). 

Even for neutron stars that emit neutrinos via modified URCA reactions in their cores, 
we predict energies and recurrence times somewhat larger than those observed. The ob- 
served superburst energies are ~ 10 42 ergs, which is roughly an order of magnitude lower 
than those in Figure 11 for modified URCA cooling. However, a substantial fraction of the 
superburst energy may be released in the form of neutrinos (Strohmayer & Brown 2002), so 
this discrepancy may not be serious. Potentially more troublesome are the relatively long 
superburst recurrence times predicted by the model. Although observations do not constrain 
recurrence times very well, the times are probably shorter than those in Figure 11. Neutron 
star models with hot cores, T corc > 5 x 10 8 K, seem to best match the observational data, 
though the M cutoff is larger (see Figure 12). Such temperatures correspond to neutrino 
cooling that is even less efficient than that from modified URCA reactions. One possibility is 
that a nonnegligible fraction of the neutron star core consists of superfluid baryonic matter. 

The above conclusions about the neutrino emission mechanism in the core assume that 
the stellar radius is the canonical 10 km. The radius of the neutron star is the only other 
parameter which we have investigated that significantly affects superburst characteristics 
at the accretion rates at which superbursts are observed. Larger stars produce superbursts 
with larger energetics and longer recurrence times (see Figure 16). Neutron stars with 
exceptionally large radii (R ~ 16.4 km) produce extremely energetic superbursts that are 
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grossly inconsistent with observations, even if the core temperature is very high. We cannot 
make such a definitive statement about neutron stars with smaller radii because we are 
unable to differentiate between superbursts from neutron stars with small radii (R ~ 6.5 
km) and moderate core temperatures (T corc k3x 10 s K) and neutron stars with moderate 
radii (R ~ 10.4 km) and high core temperatures (T core « 8 x 10 8 K). 



6. Comparison with Previous Theoretical Work 

Previous theoretical investigations of the superburst phenomenon (Brown & Bildsten 
1998, Cumming & Bildsten 2001, Strohmayer & Brown 2002, Brown 2004) generally agree 
quite well with our results. To determine the physical conditions under which superbursts 
occur, the authors use approximate ignition criteria evaluated at the base of the accreted 
layer. In this section, we compare and contrast the results of our rigorous global linear 
stability analysis with those obtained from these approximate one- zone ignition criteria. 

According to Brown & Bildsten (1998), an instability ensues when the carbon nuclear 
energy generation rate ec at the base of the accreted layer satisfies the criterion 

dec decooi { , 

dT > dT ' 

where e coo i = pKT '/£ 2 is an approximation to the global cooling rate and K is the thermal 
conductivity evaluated at the base of the accreted layer. Cumming & Bildsten (2001) and 
Brown (2004) set dlnec/dlnT = 26 and dlne coo i/dlnT = 2 and express the ignition cri- 
terion as ec > (2/26)e coo i- We presume that Brown & Bildsten (1998) and Strohmayer & 
Brown (2002) calculate the derivatives numerically. Unlike Brown & Bildsten (1998) and 
Strohmayer & Brown (2002), Cumming & Bildsten (2001) and Brown (2004) do not explic- 
itly state that they track the evolution of the carbon mass fraction Zcno when they solve 
for the equilibrium configuration of the layer. To our knowledge, these authors keep .Zcno 
constant throughout the layer, but we cannot state this with certainty. 

To conduct an accurate comparison between our global analysis and the various one- 
zone approximations, we solve for the equilibrium configuration of the accreted layer in an 
identical fashion for both methods. This ensures that any differences in the results are due 
only to the stability calculation. We find that dlne coo i/dlnT ps 2 in all of our calculations. 
Therefore, setting dlne coo i/dlnT = 2 is appropriate for all scenarios in our opinion. In 
contrast, we find that setting dlnec/dlnT = 26 is appropriate only if the core temperature 
is high (> 10 8 K) or the accretion rate is near the Eddington limit. At sufficiently high 
densities and low temperatures, dlnec/dlnT <C 26, often by many orders of magnitude, so 
this approximation is inappropriate in these situations. 
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For a given equilibrium configuration, we calculate the one-zone approximation in four 
different ways. When calculating dec/dT at the base of the accreted layer we (i) use the 
value of Zcno and d In ec/d In T derived from our equilibrium configuration, (ii) artificially set 
■^cno constant and use the value of din ec/d In T derived from our equilibrium configuration, 
(iii) use the value of Zcno derived from our equilibrium configuration and set d In ec/d In T = 
26, or (iv) artificially set Zqtsso constant and set din ec/d In T = 26. Artificially setting Zcno 
constant results in superbursts at all accretion rates for each calculation we performed. This 
is clearly incorrect since all our global stability analysis models indicate a minimum M below 
which superbursts are absent. Therefore, the criteria (ii) and (iv) above are very inaccurate 
and should be avoided. Regarding the other two criteria, we find that for neutron stars 
with high core temperatures, such as those with cores that radiate neutrinos via modified 
URCA reactions, the results from the calculations in which din ec/d In T is set to 26 and 
in which din ec/d In T is calculated self-consistently are almost identical. Therefore, the 
approximation din ec/d In T = 26 is valid in this situation. However, for neutron stars with 
low core temperatures, such as those with cores that radiate neutrinos via pionic reactions, 
the results from the two calculations differ significantly, so in this case din ec/d In T must 
be calculated self-consistently. 

See Figures 17 and 18 for a comparison between the results from our global linear sta- 
bility analysis and the one- zone approximations (i) and (iii). For the calculations in which 
both Z C no and din ec/d In T are calculated self-consistently, the results of the one-zone ap- 
proximation and the global linear stability analysis generally agree quite well. The errors 
produced by the one-zone analysis are probably smaller than those due to uncertainties in 
other parameters, such as the accretion rate, impurity concentration, elemental composi- 
tion of the crust, etc. The only significant discrepancy we find between the two models is 
the superburst recurrence times in systems with accretion rates near the critical rate below 
which superbursts do not occur. In these systems, we find from our global linear stability 
analysis that the carbon layer is marginally unstable. A large fraction of the carbon fuel 
burns stably before the full instability is triggered, resulting in a delayed burst (see Paper I), 
whereas the one- zone calculation has no way of modeling delayed bursts. These "delayed" 
superbursts have significantly larger recurrence times than those derived using the one- zone 
approximation. However, they occur over a fairly narrow range of accretion rates, so this 
discrepancy may not be serious. From our comparison, we draw the following conclusions 
regarding the application of the one-zone approximation to superburst calculations: The 
one-zone approximation is appropriate for most superburst calculations if one uses the val- 
ues of Zqno an d din ec/d In T derived from the equilibrium configuration calculation. The 
resulting superburst energies should be sufficiently accurate for all accretion rates. The su- 
perburst recurrence times should be sufficiently accurate at high accretion rates, but they will 
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significantly underestimate the true recurrence times at accretion rates near the critical rate. 
All simpler approximations, such as (ii), (Hi), and (iv) above, make serious errors for some 
choices of parameters. 

7. Summary 

In this investigation, we have carried out the first self-consistent global linear stability 
analysis of the carbon fuel on accreting neutron stars and determined the physical conditions 
under which superbursts occur. Our model reproduces the general observed superburst 
features, including burst energies, recurrence times, and the range of accretion rates at 
which superbursts occur. In contrast to normal Type I X-ray bursts, the observational 
characteristics of superbursts are very sensitive to the thermal profile of the entire crust 
of the neutron star. Consequently, superbursts can be useful probes to study neutron star 
interiors. Because our theoretical model evaluates the thermal profile self-consistently, we 
are able to explore a wide range of neutron star parameters and study the effects of each on 
the resulting superburst energies and recurrence times. 

By comparing our results with observations, we are able to set constraints on various 
neutron star parameters. We find that accreting neutron stars with highly efficient neutrino 
emission in their cores produce extremely energetic superbursts which are inconsistent with 
observations. Such neutron stars also do not have superbursts in the range of accretion rates 
at which superbursts are observed unless the crust is very impure. Stars with less efficient 
neutrino emission produce bursts that agree better with observations, while stars with highly 
suppressed neutrino emission in their cores, e.g., because of superfluidity, produce bursts that 
agree best with observations. 

If the neutron star crust is disordered, the thermal profile of the crust, in particular the 
region in which superbursts are triggered, is insensitive to the core temperature. Therefore, 
superburst energetics and recurrence times from neutron stars with highly efficient core neu- 
trino cooling mechanisms may still be consistent with observations if the crust is completely 
disordered. However, such neutron stars will have superbursts at accretion rates lower than 
the observed cutoff rate. The cutoff at M ~ O.lMEdd thus requires some other effect, e.g. 
delayed mixed bursts. 

Neutron stars with large radii (R ~ 16 km) produce very energetic superbursts that 
are inconsistent with observations, even if the core neutrino emission mechanism is highly 
inefficient. Better constraints on neutron star parameters could be made with improvements 
to the theoretical model and with more observational data. 
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All systems in which superbursts are observed and that accrete predominantly hydrogen 
have exceptionally high a values (see §5). Delayed mixed bursts, which are normal bursts 
that burn a significant amount of hydrogen and helium stably to carbon before the instability 
is triggered (Paper I), occur in systems for which 0.1 < M/M^d ^5 0-3 and produce a values 
consistent with those observed in these systems. Since the delayed mixed burst instability 
occurs above a layer of carbon that has been already produced by stable burning, it is possible 
that they leave behind sufficient amounts of unburnt carbon to fuel superbursts. For systems 
in which delayed mixed bursts do not occur, it is likely that the rp-process burns nearly all 
the available fuel to elements heavier than carbon, leaving insufficient carbon for superbursts. 
We thus speculate that delayed mixed bursts are a necessary prerequisite for the occurrence 
of superbursts in systems that accrete predominantly hydrogen. More work is needed to 
verify this suggestion. 

Previous theoretical investigations have used approximate one-zone ignition criteria to 
determine the conditions under which superburst occur. We find that the results of these 
approximate criteria generally agree well with our global stability analysis if the parameters 
used in the one-zone criteria are consistent with those derived in the equilibrium configura- 
tion calculation. Specifically, calculations in which the carbon mass fraction is kept constant 
do not reproduce the absence of superbursts at lower accretion rates, and calculations which 
set dlnec/dlnT = 26 produce inaccurate results for neutron stars with highly efficient neu- 
trino emission mechanisms in their cores. For calculations in which these quantities are 
consistent with those derived in the equilibrium configuration calculation, the only signif- 
icant discrepancy we have found is for systems that accrete at a rate close to the critical 
accretion rate below which superbursts do not occur. At these accretion rates, a substantial 
amount of carbon burns stably before a superburst is triggered. Consequently, the one-zone 
approximation drastically underestimates the superburst recurrence times. 

It is our pleasure to thank Edward Brown, Jean in't Zand, and Don Lamb for helpful 
discussions and the referee for insightful comments and suggestions. This work was supported 
by NASA grant NNG04GL38G. 
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Fig. 1. — Carbon energy generation rates as a function of rest mass density, plotted for four 
different temperatures. The composition assumed is 30% 12 C and 70% 56 Fe by mass. Note 
that at high densities, when the burning is strictly pycnonuclear, the energy generation rate 
is essentially independent of temperature. 
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Fig. 2. — Thermal and net energy flux profiles for two neutron stars with different mass 
accretion rates just before a superb urst is triggered. The parameter Z acc is the accretion 
rate normalized to the Eddington limit. The rate of accretion of matter onto the surface 
of a neutron star significantly affects the thermal state of the entire crust, as well as the 
stellar core. For this plot we assume fiducial parameter values M = 1.4 M , R = 10.4 km, 
Cf = 0.3, X out = 0.7, F out = 0.28, and Z CN0 , ut = 0.016. Additionally, we assume that the 
heavy element in the accreted layer is 56 Fe and that the core emits neutrinos via modified 
URCA reactions. These assumptions hold for all other Figures unless notified otherwise. 
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Fig. 3. — Thermal profiles for two neutron stars with different accreted gas compositions 
just before a superburst is triggered. The normalized accretion rate / acc = 0.3 for this and all 
subsequent thermal profiles. The mass fractions of the accreted gas are X = 0.7, Y = 0.28, 
Z = 0.004, and Z CNO = 0.016 for "H Accretor," and X — 0.1, Y — 0.88, Z = 0.004, and 
-Zcno — 0.016 for "He Accretor." Hydrogen burning releases much more energy per gram of 
fuel than helium burning. Consequently, the temperature of the superburst ignition region, 
10 11 gem" 2 < S < 10 14 gem -2 , is greater for hydrogen accretors than for helium accretors. 
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Fig. 4. — Thermal profiles for two neutron stars with different heavy element compositions 
in the accreted layer resulting from stable and unstable hydrogen and helium burning just 
before a superburst is triggered. "Ru" signifies a composition consisting of 30% 12 C and 



70% 104 Ru by mass, and "Fe" signifies a composition consisting of 30% 12 C and 70% 56 Fe 
by mass. The heavy element composition has almost no effect on the thermal profile of the 
neutron star crust. 
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Fig. 5. — Thermal and net energy flux profiles of the crust for neutron stars with different 
core temperatures just before a superburst is triggered. Note that a positive value of F 
denotes an outward flux. The crosses indicate the column depth at which a superburst is 
triggered. The core temperature significantly affects the thermal profile of the superburst 
ignition region, 10 11 gem" 2 < X < 10 14 gem" 2 . The flux near this region is sensitive to 
several parameters including the core temperature, and it can be either positive or negative. 
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Fig. 6. — Thermal profiles for accreting neutron stars with different values of the impurity 
parameter in the substrate, below the accreted layer. We plot three thermal profiles for a 
neutron star with a core that emits neutrinos via modified URCA reactions (curves with 
logT ~ 8.5 at logS ~ 18) and three profiles for a star with a core that emits neutrinos via 
pionic reactions (curves with logT ~ 7.4 at logS ~ 18). For each star with a given core 
cooling mechanism, the three profiles correspond to impurity parameter values Q = 100, 
5.2, and from hottest to coldest. Note that impurity scattering has a significant effect on 
the thermal profile of the crust, especially in neutron stars with cold cores for which phonon 
scattering is weak. 



-37- 




2 4 6 8 10 12 14 16 18 



logp (g/cm 2 )] 

Fig. 7. — Thermal profiles for accreting neutron stars with disordered crusts. "mURCA" 
refers to a neutron star with a nonsuperfluid core that emits neutrinos via modified URCA 
reactions, and "Pion" refers to a neutron star with a core that emits neutrinos via pionic 
reactions. The low thermal conductivity of the disordered crust makes the crust rather 
insensitive to the core temperature. Contrast these results with those shown in Fig. 5. 
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Fig. 8. — Superburst energies and recurrence times for mixed hydrogen/helium accretors 
and for helium accretors as a function of accretion rate. The core is assumed to consist of 
normal (nonsuperfluid) matter and to emit neutrinos via modified URCA reactions. The 
composition of the burned hydrogen and helium is taken to be 30% 12 C and 70% 56 Fe by 
mass (i.e. Cf = 0.3). We assume the fiducial value Cf = 0.3 for all calculations unless notified 
otherwise. 
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Fig. 9. — Superburst energies and recurrence times as a function of accretion rate for different 
values of Cf. The core is assumed to emit neutrinos via modified URCA reactions. 
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Fig. 10. — Superburst energies and recurrence times for different heavy element compositions 
in the neutron star crust as a function of accretion rate. The core is assumed to emit neutrinos 
via modified URCA reactions. 
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Fig. 11. — Superburst energies and recurrence times for two core neutrino cooling mecha- 
nisms as a function of accretion rate. The modified URCA and pion cooling models both 
assume a nonsuperfluid core. 
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Fig. 12. — Superburst energies and recurrence times for different fixed temperatures at the 
crust-core interface as a function of accretion rate. 
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Fig. 13. — Superburst energies and recurrence times for neutron stars with different impurity 
parameter values in the substrate, below the accreted layer. The core is assumed to emit 
neutrinos via modified URCA reactions. 
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Fig. 14. — Superburst energies and recurrence times for neutron stars with different impurity 
parameter values in the substrate, below the accreted layer. The core is assumed to emit 
neutrinos via pionic reactions. 
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Fig. 15. — Superburst energies and recurrence times for neutron stars with completely dis- 
ordered crusts as a function of accretion rate. 
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Fig. 16. — Superburst energies and recurrence times for different stellar radii as a function 
of accretion rate. The core is assumed to emit neutrinos via modified URCA reactions. 
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Fig. 17. — A comparison of superburst energies and recurrence times between our global 
linear stability analysis and the local one-zone approximation. "Global" signifies the results 
using our global linear stability analysis and "1-zone" signifies the results using the one- 
zone approximation for which we set dlnec/dlnT = 26. The results for which we calculate 
dlnec/dlnT self-consistently are nearly identical. The core is assumed to emit neutrinos 
via modified URCA reactions. 
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Fig. 18. — A comparison of superb urst energies and recurrence times between our global 
linear stability analysis and the local one-zone approximation. "Global" signifies the results 
using our global linear stability analysis, "1-zone A" signifies the results using the one- 
zone approximation for which we calculate dlnec/dlnT self-consistently (criterion i in §6), 
and "1-zone B" signifies the results using the one-zone approximation for which we set 
dlnec/dlnT = 26 (criterion iii). The core is assumed to emit neutrinos via pionic reactions. 



